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The entanglement entropy for the ground state of a XY spin chain is related to the 
corner transfer matrices of the triangular Ising lattice and expressed in closed form. 



In the study of the properties of quantum states the entanglement entropy plays an important role. 
It is calculated from the reduced density matrix for part of the system and measures how many states 
are connected across the interface between both parts. For critical one-dimensional systems, it diverges 
logarithmically with the size of the subsystem, and the prefactor of the logarithm involves the conformal 
anomaly c. [1,2]. This has been checked numerically for a number of solvable spin chains recently [3,4]. 
Analytical calculations exist for for the (critical) XX chain (equivalent to free fcrmions hopping on a 
lattice) [5] and for the (non-critical) XY chain in a field [6], as well as for the ferromagnetic Heisenberg 
chain [7]. In the first two cases, the limit of a large subsystem was treated by using asymptotic prop- 
erties of the Toeplitz matrices which determine the density matrix eigenvalues. For the XY chain, the 
entropy was thereby obtained as an integral over theta functions. In the present note we want to point 
out that the XY chain can also be treated by connecting it to a triangular Ising model and its corner 
transfer matrices. Moreover, a much simpler formula for the entropy can be given (at least for part of 
the parameter space) which involves only elliptic integrals. 

Our aim is to study the ground state |0 > of the XY Hamiltonian 

#=-£[(1 + 7)«+i + (1 " + 2ha*} (1) 
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where the cr" are Pauli matrices and < 7 < 1, h > 0. This operator can be obtained from 

6 = - £ [ (1 + 7K - (1 - l)<-i« +1 + 2hr^ +1 ] (2) 
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where the r° are also Pauli matrices by the dual transformation = On-i^ni t*t* +1 = cr*. The 
ground state of H will be denoted by |0 >. It is important to make this distinction because the dual 
transformation affects the entanglement properties. The bond variable a z only registers the relative ori- 
entation of the neighbouring r z -spins. 

Now it is known that H commutes with the symmetrized transfer matrix T of a triangular Ising model 
shown in Fig. 1 [8-10]. Moreover, |0 > corresponds to the maximal eigenvalue of T. The relation of the 
three coupling constants Ki — (3Ji of the triangular lattice and the chain parameters is 

7 = 1/C 3j h = (S 1 S 2 C 3 + C 1 C 2 S 3 )/C 3 (3) 

where Si = sink 2Ki and Cj = cosh 2Ki. The ordered region of the Ising system, T < T c , corresponds 
to h > 1, the disordered one to h < 1. For K$ = 0, one is dealing with a square lattice. Then 7 = 1 and 
H becomes the Hamiltonian of an Ising chain in a transverse field. This system was studied in [11] by 
noting that the reduced density matrix p for a half-chain is the partition function of a long Ising strip 
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with a perpendicular cut. If the chain length resp. the strip width are much larger than the correlation 
length, this partition function can be expressed in terms of infinite corner transfer matrices (CTM's). 



The same can be done for the triangular lattice. In this case one multiplies alternatingly CTM's with 
a central spin (Ai) and without one (Hi) as shown in Fig. 2. These quantities were studied in [12,13]. 
The result is that p has the diagonal form 

p = Kexp{ r ^ijc]c j ) (4) 

j>0 

with fermionic operators c] , cj and single particle eigenvalues 

f (2j + 1) s, T>T C 

ij = (5) 
1 2js, T<T C 

The factor e which sets the scale is e = tt I(k')/I(k) where I(k) denotes the complete elliptic integral 
of the first kind and kl = Vl — k 2 . The modulus k, finally, is the Onsager parameter for the triangular 
lattice and given by eqn. (6.4.16) in Ref. [14]. In terms of the chain parameters one finds, using (3) and 
choosing K\ = K 2l that 



Jh 2 + 7 2 - 1 / 7, h<l 

k = { / (6) 

7 / ^/h 2 +j 2 -l, h > 1 

To obtain the reduced density matrix p for the XY chain, one notes that the dual transformation inter- 
changes order and disorder. Therefore one has to interchange the expressions for the eigenvalues in (5) 
to obtain 

f 2j e, h< 1 
{ (2j + 1) e, h>l 

In the special case of the transverse Ising chain (7 = 1) this spectrum coincides with the one obtained 
directly from (2) and (4). The eigenvalue e = describes the (asymptotic) degeneracy of the density 
matrix cigenfunctions for h < 1 connected with the spin inversion a x — > —a x . 

This is valid for purely ferromagnetic interactions as well as for antiferromagnetic competing interac- 
tions in the Ising model. If we confine ourselves to the first case, h and 7 satisfy h 2 +-f 2 > 1 and k is real, 
< k < 1. Thus one never crosses the "disorder line" h 2 + 7 2 = 1 in the 7 — /i-plane. In the terminology 
of [6] one is in the regions (la) and (2). In these regions we have completely determined the reduced 
density matrix. Comparing with Ref. [6] one sees that the eigenvalues ej correspond to the quantities 
there. Moreover, e is precisely the quantity ttto, as can be seen by evaluating the integrals for t . The 
two approaches therefore give the same density-matrix spectrum, as they should, but the way via corner 
transfer matrices provides more physical insight. 

Given the thermal form of the reduced density matrix, the entanglement entropy S = —tr(p Inp) can 
immediately be written down 

S = ^ In [1 + expi-ej)} + ]T 1 = InZ + U (8) 

j j v\ 3) 
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This expression has been evaluated numerically and analysed near criticality (where e goes to zero) in 
[2]. However, one can also obtain a closed form by introducing the elliptic nome 



q = exp(-s) = exp(—n I(k')/I(k)) 



(9) 



and invoking elliptic function identities. In the case ft > 1, one finds the following formula for Z from 
equ.(16.37.3) in [15] 



11(1+/ 

3=0 



2j+l 



\ _ I Ife? \l/24 

' { k 2 k'2' 



while for the sum in the "internal energy" U, equation (16.23.11) in [15] leads to 



3=0 



1 



+ g 2 J'+! 24 



[1 - (1 -2fc 2 )(2/(fc)/^) 2 



Similar formulae hold for the case ft < 1. In this way one arrives at the results 



S =T2 



_ k 2 AI(k)I(k') 



+ ln2, h < 1 



(10) 



(11) 



(12) 



*=24 



ft > 1 



(13) 



where the In 2 in (12) arises from the eigenvalue £o = 0. Near criticality (ft — ► 1, k — > 1) these formulae 
give the logarithmic divergence 



(14) 



found in [2] where the prefactor corresponds to the value c = 1/2 of the central charge. For k — > 0, the 
brackets vanish in both cases. If ft < 1, this happens as one approaches the disorder line. Then |0 > 
becomes the superposition of two product states [10] and S goes to In 2. The corresponding divergence 
of the sj for j =/= 0, which leads to a collapse of the spectrum of p, has been seen in numerical calculations 
for small systems [16]. A particular case is 7 = 1, where one approaches for /nOa superposition of the 
two states with all spins in the positive resp. negative x— direction. For ft > 1, the parameter k vanishes 
as ft — > 00. Then |0 > becomes the simple product state | + + + > and the entropy vanishes. 



All these results refer to a chain with open ends which is divided into two parts. These two parts have 
one point of contact. If one cuts a segment out of a long chain or a ring, which is the situation studied in 
[6], one creates two such interfaces. Since the density matrix cigenfunctions for the low Sj are localized 
near an interface, one obtains a double degeneracy of the Sj with j =/= if the segment is longer than the 
correlation length. This degeneracy was also noted in [6] and gives a factor of two for the brackets in 
(12) and (13), while the logarithm is unchanged. 

We have not discussed the region j 2 + h 2 < 1 inside the disorder line (the region (lb) in [6]). The 
CTM's for this case have not been studied so far. However, on the line h = one can relate the Hamilto- 
nian of the XY chain to the transfer matrix of a doubled square lattice, i. e. an eight-vertex model at the 
decoupling point [17,14]. For this case the CTM's are known again [14]. One then obtains Sj = j e and 
the value for k is k — (1 — 7)7(1+7). This corresponds again to the result for 7rr in [6]. Since now even 
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and odd multiples of e enter, the entropy for a chain cut in half becomes the sum of the two expressions 
(12) and (13). The divergence of S in the isotropic limit 7 — > 0, k — > 1 is therefore also stronger than 
(14) by a factor of two, in accordance with the value c = 1 for this critical point. 
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FIG. 1. Portion of a triangular lattice forming the row transfer matrix T 




FIG. 2. Construction of the reduced density matrix p from six triangular corner transfer matrices. Along the 
heavy line the lattice is cut open. The notation follows Ref. [13]. 
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